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ABSTRACT 

We analyse new spatially resolved integral field spectroscopic H- and K-band data at a resolu- 
tion of 0.3" (100 pc) and re-analyse interferometric CO(2-l) line observations of the prototypical 
merging system Arp 220. We find that the majority of the K-band luminosity is due to a lOMyr 
old starburst, with a significant contribution from an underlying > 1 Gyr ol d stellar population, 



and a small contribution from stars < 8 Myr old. The ICalzetti et al.l ([20001 ) reddening law pro- 
vides the best fit to photometric datapoints spanning 0.45 /im - 2.12 fim. Furthermore, estimates 
of the bolometric luminosity from IRAS fluxes in conjunction with our stellar population analysis 
indicate that we observe less than 10% of the emitted K-band light. The stellar and CO(2-l) 
kinematic centre of the western nucleus coincides with the compact hot dust emission, indicat- 
ing that the latter marks the centre of the gravitational potential. In the eastern nucleus, the 
CO(2-l) data are well matched by a model in which the gas orbits around the peak of the dust 
emission. This, and the similarity of the K-band tracer kinematics, shows that despite the irregu- 
lar morphology, the eastern nucleus is also a kinematically coherent structure. Comparison of the 
extinction map and EWco and EW_b,- 7 maps indicates that the lower half of the eastern nucleus 
is significantly more extincted than the upper half, suggesting that the lower half is buried in the 
larger-scale gas disk. 

Subject headings: galaxies: active — galaxies: individual (Arp 220) — galaxies: interactions — 
galaxies: kinematics and dynamics — galaxies: starburst — infrared: galaxies 



Introduction 



Local Ultraluminous Infrared Galaxies (ULIRGs, with Ltr > 10 12 L m , Sanders fc Mirabel 19961) have 



received considerable atten tion since their discovery by IRAS ( Neugebau er et al. 1984 ) . Th e y are almos t 



exclusively major mergers (panders et al.lll988t ISanders fc Mirabel 



1996 



Veilleux et al 



2002; [JogeeJ [2004), 



and widely h eld to constitute a key epoch in the transformation of gas-rich spiral galaxies to elliptical 



galaxies (e.g. Hopkins et al. 20061 ). Although not making a significant contribution to the cosmic energy 



budget at current times, th eir relative importance increases towards higher redshifts (e.g. ICaputi et al.ll2007 



Magnelli et al 



2009 



201(3). Thus the study of local ULIRGs is important for several reasons; to properly 



understand the major merger process believed to trigger qua sar activity and res ult in elliptical galaxies; and 
as local analogues of the most luminous high-z galaxies fe.g. lEngel et al.ll2010bl ). 

Arp 220 (L IR = 1.4 X 10 12 L g . ISoifer et all 19871 ) is the closest ULIRG representative (D = 73 Mpc; z = 0.018, 
l"=352pc) and thus a unique laboratory for the investigation of this important class of galaxies. As such, 
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it has been extensively studied and is widely regarded as the prototypical ULIR G. It is an advanced merger; 
tidal tails and distortions a re observed at optical wavelengths and in HI emission |Arp|ll966t iJoseph fc Wrightl 



1985; Hibbard et al 



tions (Armus et al 



199a iHeckman et all 119961 ; iMcDowell et al 



2000) . There is evidence f o r a galactic-scale outfl ow from Ha and soft X-ray observa- 



1997; iDownes fc Solomon! Il998t ISakamoto et al.l Il999t iDownes fc Eckarti 120071) and radio ( Mundcll et al 



2003). In the centre, two nuclei are dis- 



cernible with a separation of ~ .4kpc in the near- infrared (Scoville et al. 19981) . millimeter ( Scoville et al 



2001 ) regimes. There also is a more extended (kpc-scale), rotati ng gas disk (jDownes fc S olomon 



1998; 



Sakamoto et~alll999ll2009llMundell et al.lbOQlUwiedner et aLll2002h . Arp 220 is very gas-rich: IScoville et al 
( 19971) measure a molecular gas mass of ~ 9 x 1O 9 M0, concentrated in the central 750 pc - implying an aston- 
ishing molecular gas surface density of 5 x 10 4 M Q pc~ 2 . It therefor e comes as no surp r ise that obscuration 



is sev ere for this system; estimates for Ay range from 50 to 1000 (ISturm et al 



1996 



Downes fc Solomon 



19981 ). and even at near-infrared wavelengths obscuring dust lanes a re visible (jScoville et al.lll998l ). Observa- 
tions of CO (3-2) and th e 860 /j,m continuum |Sakamoto et al. 2008) and CO (6-5) and the 435 /mi continuum 
(jMatsushita et al.l 120091 ) indicate substantial dust optical depths even at submillimetre wavelengths. These 
extreme levels of extinction have hindered a definitive determination of t he power source of Arp 220; a signif- 
icant contribution to its prodigious luminosity comes from a starburst ( Sturm et al. 19961 ; Lutz et al. 19961 
1998tlGenzel et al.lll998l ). but there is also sign ificant (and possib ly dominant) dust emission from heavily ob- 
scured central star formation or perhaps AGN (jSpoon et al ] l2004h . The presence of a deeply dust-ensh r ouded 
major nuclear power source is confirmed by millimetre observations; most recently Sakamoto et al. ( 20081) 
estimate the bolometric luminosity of the 50-80 pc core within the western nucleus to be at least 2 x 10 11 Lq 
and possibly as much as ~ 10 12 L Q . The hig h obscuration has prevented confirmation or reject ion of a possible 



active nucleus through X-ray observations (Iwasawa et al 
colum n density is of order Njj 



2001 



10 25 cm 2 (Iwasa wa et al 



2001 



2005 



Clements et al 



Saka moto et al.l 12008). 



2002) ; the obscuring 



Downes fc Eckart 



([2007Q interpret the very compact dust continuum emiss ion in the western nucleus as being heated by a 
black hole accretion disk, whereas Sakamoto et al. ( 20081) contend that the starbursting west nuclear disk 
must have in its center a dust enshrouded AGN or a very young starburst equivalent to hundreds of super 
star clusters. 

Since Arp 220 is the ULIRG most frequently used as a starburst gala xy template, it is vita l to understand 



its power source and star formation history. A photometric analysis by IWilson et al.l (|2006l ) of HST UB VI 
observations showed that the central cluster population divides into two groups; one with ages < 10 Myr, 
and one with ages ~ 300 Myr. However, it must be kept in mi nd that due to photometric u ncertainties, 



the age/reddening degeneracy limits the accuracy of such work. iRodrfguez Zaurin et al.l (|2008l ) use optical 
spectroscopic observations of the extended diffuse light along three slit positions for a stellar population 
analysis; they find the optical spectrum to be dominated by an intermediate-age stellar population (0.5- 
0.9 Gyr), and a young stellar population (< 100 Myr) which has an increasing contribution towards the 
central part of the galaxy. Besides the extinction limitations afflicting all optical- wavelength observations of 
dusty galaxies, it must be cautioned that only one of their slits covers the central region, from which the vast 



Parra et al. ( 


2007) (see also 


Rovilos et al. 


2005; 



2006) investigate radio supernovae and supernova remnants in the nuclear region of Arp 22, finding that they 
are indicative either of a radically different stellar initial mass function, or a very short, intense burst of star 
formation ~ 3 Myr ago. 

Here, we use adaptive optics near-infrared integral field spectroscopy data and interferometric mm CO(2- 
1) and continuum observations, to investigate the star formation history, the extinction, and the nature and 
orientation of the two nuclei. 
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In 12 we introduce the data, and data analysis procedures. In Sj3] we investigate the star formation 
history in the SINFONI field of view, and in S0]look at the extinction in the near-infrared regime. We then 
focus on the western (jj5]) and eastern (^BJ) nucleus. We put this together in a coherent picture of the nuclear 
region ($7]), before summarising and concluding in ^Hl 



2. Observations and Data Processing 
2.1. SINFONI Data 

2.1.1. Observations and Reduction 



Observations of Arp 220 were performed on the nights of March 7th and April 18th to 21st 2007 on Cerro 
Paranal, Chil e, at the Very Large Tel escope (VLT) with SINFONI. SINFONI is a near-infrare d integral field 
spect rometer ( Eisenhauer et al. 20031) which includes a curvat ure-based adaptive optics system ( Bonnet et al 



20041 ) and can operate with the VLT's laser guide star facility ( Bonacini Calia et al J2006t iRabien et al . 2004) 



The laser guide star was used for these observations, without tip-tilt correction due to lack of a suitable tip- 
tilt star. SINFONFs 0.05" x 0.10" pixel scale was used, giving a 3.2" x 3.2" field of view. The science data and 
standard star observations were reduced using the SPRED software package (jAbuter et al.l 120061 ). Telluric 
correction and flux calibration of the reconstructed data cubes were performed using these standard star 
frames. Fig. Q] shows spectra from three different positions of the Arp 220 nuclear region. 



2.1.2. PS F Estimation 



Several avenues exist to estimate t he poin t spread function (PSF hereafter) of adaptive optics data. 
Here we employ the method outlined by iDaviesI (|2008al) and used e.g. bv lMueller Sanchez et al.l ([20061 ) and 
Engel et al.l (|2010al ). which proceeds by comparing the data for which the PSF is not known, denoted Ii ow 
here, to higher- resolution data with known PSF (Ihigh)- Since an observed image is the convolution of the 
instrumental PSF with the intrinsic on-sky image, I b s — lintr ® PSF, the PSF of the AO data can be 
estimated by finding the broadening function B which, applied to the higher-resolution image, yields the 
best match to the AO image; Ii ow = Ihigh ® B. Since I intr is the same for both Ii ow and Ihigh, the PSF of 
how can be estimated by convolving the PSF of the high resolution image with the broadening function B; 
PSF low = PSF high ® B. 

In our case, we use HST NICMOS 2.16/im observations retrieved from the HST archives. Fig.[2]shows 
our SINFONI data, the HST data (rotated and resampled to the SINFONI pixel scale), and the broadened 
HST image. We find the PSF to be well represented by an almost perfectly symmetric Gaussian with major 
axis oriented 4° east of north and F WHM 0.30" x 0. 31"- a good performance for laser guide star assisted 
adaptive optics without a tip-tilt star |Daviesll20o"8bT ). 



2.1.3. Spatial binning 



In order to achieve a mininum signal-to-noise , we spatially binned e ach spectral plane of our reduced 
data cube using an optimal Voronoi tessellation ([Cappellari et al.l 120031 ). which bins pixels together into 
groups by accreting new pixels to each group until a pre-set signal-to-noise-cutoff is reached. 
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The signal-to-noise cutoff was chosen such as not to compromise spatial resolution in the central regions 
(i.e. such that no binning was performed there), but at the same time extending the region in which mean- 
ingful analyses can be performed. As a result of this, the outer regions in e.g. the velocity maps appear in 
blocks, rather than individual pixels. 



2.1. 4- Extracting Stellar Kinematics 



The K-band covers the stellar CO absorption bandheads longward of 2.29 /im, whose sharp blue edges 
are very sensitive to stellar motions. Here we utilise CO 2-0 and CO 3-1 to derive 2D maps of the stellar 
velocity and dispersion. After normalising the spectra with respect to a linear fit to the line-free continuum, 
we subtract unity and thus set the continuum level to zer o. A stellar templat e spectrum, of the M3III star 
HD 176617, was prepared analogously. As we discuss in lEngel et al.l (|2010al ). with this method - setting 
the continuum level to zero - the choice of template star is not critical, with even significantly mis-matched 
templates resulting in errors of at most only a few percent. Here, the template is matching the bandheads 
quite well. We furthermore note that finding a single template (or a composite template) that provides an 
excellent match everywhere is almost imp ossible to find, due to the large spatial variations in EWco seen 
in this system and also e.g. in NGC6240 ([Engel et al.ll2010al ). 



In order to obtain the 2D velocity and dispersion maps, we fit the template to the spectrum of each 
spatial pixel (spaxel hereafter) by convolving it with a Gaussian, varying the Gaussian parameters until a 
best match is achieved. The Gaussian position and width then yield the stellar velocity and dispersion at 
that spaxel. During the fitting, the spectral regions covering the steep edges of the bandhead at all expected 
stellar velocities are given five-fold weight, in order to focus the fit on the kinematics rather than details 
in the spectrum. The spectral ranges u s ed to compute the CO2-0 bandhead equivalent width (EWco 
hereafter) are those of iForster Schreiberl (|2000|); these ranges are also used by the stellar synthesis code 
STARS employed later. We note that we do not attempt to extract the Gauss-Hermite terms /13 and /14, 
since the signal-to-noise across most of our spectra is not sufficient to include these in the fit. 



The resulting final maps of stellar velocity, dispersion, and EWco are displayed in Fig. |3l 



2.1.5. H- & K-band Gas Tracers 

We derive emission, equivalent width, velocity, and velocity dispersion maps of H 2 1-0S(1) (A = 2.1218 /im), 
Hel (A = 2.0581 /mi), Paa (A = 1.8751 /mi), [Fell] (A = 1.6440 /mi), Br 7 (A = 2.1662 /mi), and Br<5 (A = 1.9445 /mi), 
by fitting Gaussians to the lineshapes. The results are shown in Fig. |U 



2.2. Plateau de Bure Interferometer Data 



We also analyse 12 CO(J = 2-l) line emission and 1.3 mm dust continumm emission observations from 
the IRAM millimetre interferometer, located on the Plateau de Bure, France, at an altitude of 2550 m 
( Guilloteau et al. 19921 ) . These observations were first analysed and published by Downes fc Eckart ( 2007 ) , 
and we refer the reader to this publication for the technical details. These authors find evidence for a hot 
dust source in the western nucleus, which they interpret as being heated by an AGN accretion disk; the best 
evidence for a black hole in the central region of Arp 220 to date. Since the absolute astrometry between 
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the PdBI and the SINFONI datasets is not precise, we anchor the SINFONI and PdBI data with respect 
to each other by assuming that the flux peaks of the western nucleus are coincident. This assumption is 
justified by the fact that in the western nucleus, all the K-band tracers coincide at the same location (unlike 
in the eastern nucleus), and hence it is plausible to expect that the cold gas emission also peaks at the same 
position. We derive a velocity map by fitting Gaussians to the line shapes, the result is displayed in Fig. [5] 



3. Star Formation History &: Stellar Populations 



In this Section, we investigate the star formation history, a nd the composit i on of the stellar pop- 
ulation, using observables and the stellar synthesis code STARS ([Sternberg 11998c ISternberg et al.l 12003 ; 
Fdrster Schreiber et al.ll2003l bavies et aDl2003l b005l bood hooii) . In the following, we use spatially aver- 
aged measurements to investigate the composition of the stellar population. Since we are summing over a 
large spatial area, we can be sure to cover a large number of individual star clusters (which may have started 
forming stars at slightly different times), and hence to be probing the average star formation properties. 
Our primary diagnostics are the EWco and the Br7 flux and EWs r7 . The EWco map is shown in the last 
panel of Fig. [3l the Br7 maps are displayed in Fig. [6] In Fig. we show the evolution of EWco, EWs r7 , 
and Lboi/Lx for a number of starbursts with different decay timescales. We furthermore use dynamical 
mass estimates of the nuclei (details can be found in ij5.3l and ij6.2|) and predictions of the mass-to-light 
ratios of different stellar populations (Fig. [7J. In what we follows, we show that a single stellar population 
cannot account for the observed diagnostics. Instead we find that three distinct populations are required; a 
10 Myr old starburst accounting for over half the observed luminosity, a significant contribution from an old 
(> 1 Gyr) stellar population, and a rather smaller population of very young stars, responsible for the Br7 
emission. 



3.1. A 10 Myr old Starburst... 

We observe EWco to be > 14 A everywhere, and > 16 A on the nuclei; as can be seen from Fig. [Jj only 
an instantaneously decaying ('delta'-) starburst 10±2Myr old can produce such large values of EWco- Do 
these stars account for the entire stellar luminosity and mass, or are there significant contributions from 
other stellar populations? We can estimate the luminosity contribution of any stars not belonging to this 
10 Myr old starburst population by noting that the theoretically possible upper limit to EWco is 17 A. This 
allows us to calculate the luminosity contribution from a stellar population with smaller EWco that may 
be diluting it. We note that this is strictly speaking only an upper limit, since we do not know the intrinsic 
EWco of the starburst population. However, this intrinsic value likely is quite close to 17 A, since this is 
what is measured in between the two nuclei, where the dilution from any putative older stellar population 
associated with the nuclei would be minimal. Hence our 'upper limit' to the non-starburst luminosity is in 
fact more likely to be a reasonable estimate. 



3.2. 



& a younger stellar population? 



Both the presence of radio supernovae ( Parra et al and our Br7 maps clearly indicate the presence 

of young stars - only very young (< 8 Myr) stars are hot enough to excite hydrogen sufficiently to emit Br7 
flux. An important question is whether this population contributes to the K-band continuum. We thus first 
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investigate the possibility that a population of stars younger than lOMyr may be diluting the CO bandheads 
of the lOMyr starburst population; these stars would have EWco=0A (Fig. [7]). In this case, the younger 
stars would contribute ~ 10% of the K-band luminosity. Below O, we calculate that if the assumption of 
younger (< lOMyr) stars diluting the CO bandheads is correct, then we are only seeing « 20% of the actually 
emitted (no reddening-correction) K-band luminosity. For a 10 Myr old instantaneously decaying starburst, 
STARS predicts a K-band mass-to-light ratio of wO^Mq/Lq (Fig. [7]). Since we know the luminosity of 
the 10 Myr old stellar population under this scenario to be ^90%, we can therefore calculate its expected 
mass (correcting the measured K-band luminosity for the western nucleus for the missing light). Using 
our dynamical mass estimate f t)5.3l and assuming a 10% gas fraction), we can then calculate the expected 
mass-to-light ratio of the non-starburst stellar population of the western nucleus; we find that it must be 
M/Lr- ~ 75M Q /L Q . The same calculation for the eastern nucleus yields an even larger M/Lr- ~ 99 Mq/L q . 
The expected M/Lr- for stars < 10 Myr old is <O.1M /Lq - a divergence of three orders of magnitude. 
We therefore reject the assumption of a younger stellar population diluting the bandheads of the 10 Myr old 
starburst population. 



3.3. 



&: an older stellar population! 



Stars older than ~ 20 Myr have an EWco ~ 12 A; a fractional luminosity contribution of 26-44% would be 
required to achieve the bandhead dilution (Table[TJhsts the results for both the western nucleus and the entire 
nuclear region). In this case, we are only seeing ss7% of the actually emitted K-band luminosity (@J this 
is different to the fraction calculated above under the assumption of a younger stellar population, since the 
theoretically expected Lboi I Lk is different in this scenario) . Calculating the expected mass-to- light ratio for 
the older stellar population in the same way as above, we find that it must have M/Lr- ss 7.4M /L - which 
is what is expected for a several Gyr old stellar population. The same is found for the eastern nucleus, here 
an M/Lj<- « 9.8 Mq/Lq is required. We therefore conclude that the assumption of a significantly (> 1 Gyr) 
older stellar population alongside the 10 Myr old starburst in the central kpcs of Arp 220 is correct. We note 
that we still require a third population of stars younger than ~7Myr (which would have EWs r7 ^600A, 
Fig-©, to account for the Br7 emission. However, their luminosity contribution would be minimal; their 
intrinsic EWg r7 ~ 600 A has to be diluted by a factor > 75 to reach the measur cd EW B <8A. 



This is in good agreement with what iRodriguez Zaurin et al.l (]2008f ) find based on optical spectroscopy 
of the central regi on of Arp 220 ; name ly near-equal contributions from a < 0.1 Gyr and a 0.5-0.9 Gyr old 
stellar population. IWilson et al.l (|2006l ) derive photometric age estimates of clusters from HST UBVI data 
and find two age classes; < 10 Myr and 70-500 Myr. 



4. Reddening Si Extinction 
4.1. Reddening 

Like most ULIRGs, Arp 220 contains significant amounts of dust. In the K-band, the effect of this is 
mostly a wavelength-dependent extinction due to scattering and absorption, with longer- wavelength radiation 
being less severely affected. The effect of this is 'reddening', a steepening of the continuum level towards 
longer wavelengths. Another effect pertaining to the presence of dust is the re-emission of dust-absorbed 
stellar light. However in the H- and K-band this would only be observable if very hot dust were present. 
Since any continuum emission due to hot dust would lead to dilution of the CO absorption bandheads, and 
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since we are seeing absorption depths very near the theoretically possible maximum of ' 
the possibility of any significant hot dust emission. 



17 A, we can exclude 



When accounting for the obscuring effects of dust, the dust is most commonly assumed to either be 
spatially uniformly mixed with the stars ('mixed model'), or to form a scr een between the obse rved stars 
and the observer. Or one can use the empirically derived reddening law of ICalzetti et al.l ([20001 ). In cases 
of strong obscuration, the mixed model may be insufficient to account for the reddening, since it saturates 
beyond a certain level of extinction. In order to determine which method is appropriate for Arp220, we 
obtained archival HST images spanning a wavelength range of 0.45 -2.12 /xm, and extracted photometric 
data points at a number of locations using 0.15" radius apertures. We then reddened a synthetic stellar 
population spectrum (corresponding to the star formation history determined in Sj3|), using each reddening 
prescription outlined above, until it fit those data points. The results can be seen in Fig. |8]- as expected, 
the mixed model was insuf ficient to account for the strong extinction of Arp220. The screen model results 
in a decent match, but t he ICalzetti et al.l (|2000l ) reddening law clearly produces the best results. We note 
that lEngel et al.l (|2010al) arrive at the same conclusion for NGC 6240. Here, the screen extinction law also 
yields a somewhat lower F t s /F em than the ICalzetti et al.l (|2000r ) reddening law. 



In o rder to quantify the r eddening in Arp 220, we fitted a stellar template to the line-free continuum, 
using the Calzetti et al. ( 2000f) reddening law. In Fig.O we show the resulting map of F i, s /F ern at 2.12 /im, 
and of the corresponding optical extinction Ay\ obtained by calculating F b s /F em at 0.55 /im and converting 
to extinction via Ay = —2.5\og(F i, s /F em ). 

We also derived the reddening using molecular hydrogen emission lines. For case B recombination 
at T — 10,000 K and n e =10 4 cm -3 , the line flux ratio of Paa and Br7 is 12.2; by c omparing this to the 
actually observed line ratio, the relative reddening can be deduced. Again using the Calzetti et al. ( 2000h 
reddening law, we thus measured F f, s /F em and derived Ay, both shown in Fig. [TO] As can be seen, the 
two extinction measurements agree well both in regard to strength and spatial variation of the extinction. 
The 1.3 mm continuum map of lDownes fc Eckartl ([20071 ) could give us valuable information about the spatial 
distribution of the obscur ing dust, however unfortu nately only two barely resolved point sources are detected 
at their sensitivity level. iDownes fc Eckartl (|2007[ ) estimate the optical depth for the western nucleus to be 
T(vjv(l-3mm) > 0.7; adopting their assumptions also for the eastern nucleus, we estimate tstv(1-3 mm) > 0.15; 
this is by extension also an upper limit of the off-nuclear r(1.3mm). For rcxA~ 2 , r^O.15 implies that all 
emission shortwards of 500 /im is optically thick. There appear to be two different regimes: the extremely 
opaque central ~50pc which are fully obscured at near infrared wavelengths (r = 1 at 500 /im), and the outer 
parts where we can detect reddened near-IR emission (and where we find r = 1 is at ~ 2.5 /im). This plausibly 
explains why even after applying the reddening correction, we do not recover all the near-IR luminosity: we 
are simply not seeing all that i s emitted in the ce ntral region. We note that this is qualitatively similar 
to what we found in NGC 6240 (jEngel et al.ll2010a|) : over most of the disk the near- infrared obscuration is 
moderate, but increases dramatically in the central tens of parsecs of each nucleus. 



4.2. Extinction 

The significant reddening quantified in §4. H is indicative of a substantial amount of stellar light emitted 
in the near-infrared in the central region of Arp 220 which is completely absorbed, and re-emitted at longer 
wavelengths. Here, we attempt to quantify the amount of 'missing light' in the K-band, by comparing the 
infrared luminosity, calculated from IRAS flux measurements, to the stellar bolometric luminosity calculated 
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from o ur measured Lk and the \-h, i/Lk expected from our stellar population analysis ($3]& Table [TJ). We 
use the lSanders et al.l (|2003) IRAS fluxes from the Revised Bright Galaxy Sample, and calculate the infrared 



flux via F Zi{ ,8-iooo^m = 1.8 x 10 



(13.48 xfi 2 



flux density in Jy at 12 /xm etc (panders fc Mirabel 119961 ) 
Lf r = 3.127 x 10 7 D 2 FfR (D in pc). This yields L /flj8 



5.16xf 25 + 2.58xf 60 + fioo)Wm- 2 , with f 12 the IRAS 
We then convert this t o infrared lum i nosity , 

"~ d201dh 



Nardini et al 



j/_r = o.iz/ x iu' u~ rjB ^e (u m pcj. mis yieius W^s-lOOO/im — 1-37 X 10 12 L Q . 
perform a 5-8 /xm spectral analysis of Arp220 using the Infrared Spectrograph onboard Spitzer; they find 
a 17% AGN contribution to the bolometric luminosity of Arp220, which we need to correct for in order 
to find the stellar bolometric luminosity. However, this is largely offse t by the translation of infrared to 
bolometric luminosity, which for ULIRGs is typically Lboi= 1.15 x Lm (jKim fc SandersHl998T ). This yields 
a total stellar bolometric luminosity of 1.31 x 10 12 L Q . We then measure the K-band luminosity of Arp220 
using the archival HST NICMOS images already utilised in 32~L2l finding L K = (1.75±0.1) x 10 9 L Q (for 
comparison, our SINFONI data yield L K = (1.34±0.07) xl0 9 L f 7 ) for the SINFONI FoV, im plying that - 80% 
of the total Lk are emitted in the nuclear region - cf. IWvnn- Williams fc Becklinlll993l) . This results in a 
Lboi/LK of ~ 750, whereas our earlier stellar population analysis (j£3]) requires ~ 150 (theoretically expected 
for 10 Myr old starburst plus younger stellar population with the derived relative luminosity contributions; 
cf. Table [TJ) or ^50 (10 Myr old starburst plus older stellar population) - implying that in the K-band we 
are missing a factor of ~ 5 or ~ 15 of light, respectively. As we show in Sj3j the assumption of a population 
younger than 10 Myr is inconsistent with our dynamical mass estimate, and we can therefore exclude this 
alternative. We thus conclude that in the K-band, we ar e only observing ~ 7% of the a ctually emitted (no 
reddening-correction) stellar light. This is consistent with iRodriguez Zaurin et al.l (|2008r ). who find that the 
bolometric luminosity derived from optical wavelength observations is an order of magnitude smaller than 
the mid- to far-infrared luminosity, and are led to conclude that 'most of the on-going star formation in the 
nu clear region is hidde n by dust'. Further support for our result is lend by the supernova rate estimates 
of lRovilos etafl fcood) : these authors derive SNR~0.7yr 1 for the western nucleus. STARS predicts a 



supernova rate of 10 10 L K 1 yr 1 for a 10 Myr old instantaneous starburst, which for our corrected K-band 
luminosity, of which ~ 75% are attributed to the 10 Myr old starburst, implies a supernova rate sa 1.0 yr -1 . 



5. Western Nucleus 

5.1. Stellar & CO(2-l) Kinematics 

As Fig. [3] shows, the stellar kinematics exhibit a clear rotational signatures around the western nucleus. 
In order to extract a rotation curve and find the kinematic centre of the stellar rotation, we fit an inclined 
disk model to the velocity map (Fig. ITT)) . The best-fitting model has an inclination q = 0.62±0.06, and a 
position axis (PA hereafter) of 9.8±0.4deg south of west, with average residuals of 12.2kms _1 . We also 
derive a dispersion profile by azimuthally averaging the 2D dispersion map. Note that we do not suggest 
that the stars in the western nucleus are moving in a thin disk - we simply parameterise the rotation field 
in terms of circular orbits with a fixed centre, PA, and inclination, in order to derive a rotation curve and 
kinematic centre. 

The CO(2-l) also displays a regular rotation pattern around the western nucleus, which we model 
analogous to the stellar velocity analysis. The best-fitting model (Fig. fT2|) indicates an inclination qw0.56 
and PA«3deg south of west, with the kinematic centre coincident with that of the stellar velocity field 
within the uncertainties. The stellar and CO(2-l) kinematics in the western nucleus thus seem to agree well. 
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5.2. Kinematic Centre &c Hot Dust Emission 



Downes fc Eckartl (|2007f) detect emission from a hot, compact dust ring just south of the CO(2-l) 



emission peak of the western nucleus, which they interpret as being heated by an AGN accretion disk. 
We locate the position of the continuum emission using the relative position of the dust emission peak 
and the CO(2-l) emission, and assuming that the peak of the stellar continuum in the western nucleus is 
coincident with that of the CO (2-1) emission (as outlined in §2.2j) . Intriguingly, the putative AGN position 
thus derived is coincident with the kinematic centre of the stellar and CO(2-l) rotation found in tj5.ll This 
strongly suggests that th e stars in the western nu cleus are moving in a gravitational potent ial with either a 



supermassive black hole ([Downes fc Eckartl 120071 ) or an extremely dense, young starburst ([Sakamoto et al 
2008) at its centre. 



5.3. Dynamical Mass Estimate 



In principle, our data are of high enough quality to derive the dynamical mass via Jeans modelling. 
However, as we found in §4] Arp220 is severely affected by extinction in the K-band. Since we do not 
know the relative three-dimensional distribution of the obscuring material, it is impossible to derive the 
intrinsic kinematics from the observed two-dimensional line of sight projections. We therefore h ave to content 



oursel ves with a simpler estimate of the dynamical m ass. We adopt the approa ch taken by iBender et al 



(|1992l ) , which is based on the ten sor virial theorem (jBinnev &: Tremaind 119870 , and refer the reader to 
appendix B of IBender et al.l (|1992T ) for a detailed discussion. According to this, the appropriate dynamical 
mass formula for the relative stellar velocities and dispersions of the western nucleus is Md yn = 1.12x3 XOq X 
R/G. Using the stella rotation and dispersion profiles derived in £ 15.11 this yields an enclosed dynamical 
mass out to - 100 pc (180 pc) of 6.3 x 10 9 A/ Q (9.0 x 10 9 M Q ). 

Since the CO data probe significantly deeper th an the K-band stell ar kinematics, we use these to obtain 
another mass estimate. We use the code described in ICresci et al.l (|2009h to model the CO emission as arising 
from a thick rotating disk, taking beam smearing effects into account. We achieve the best match to the 
observed flux, velocity, and dispersion maps with a mass distribution in which the majority (> 80%) of the 
mass is concentrated within a radius of ~ 0.15" (50 pc). An even more concentrated mass distribution (point 
mass) also gives a good match (although the more extended distribution is preferred) - the 0.3" beam size 
prevents a definitive differentiation here. The total mass predicted within a radius of 100 pc is 1.6 x 10 9 Lq, 
although we note that our model underpredicts the central dispersion by ^60kms, and hence this value 
probably is somewhat too low. Since our modelling cannot definitively exclude a significant mass distribution 
from a central point mass (i.e., black hole), and since the presence of such a non-luminous mass would have 
ra mifications for our ana lysis in ^3[ we note here that the black hole mass expected from the Ms^-cr-relation 



of iTremaine et al.l (|2002n is 1.4 x 10 8 M Q ; i.e. less than 10% of the total dynamical mass. 



Downes fc Solomon! (|1998l ) estimate a gas mass for the western nucleus of 6xlO 8 M , based on 1.3 mm 



dust emission, implyin g a gas fraction of ~10%, which agrees well with what is expected for local ULIRGs 
and starburst galaxies ([Hicks et al.ll2009l and references therein). 
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6. Eastern Nucleus 

6.1. Stellar &: Gas Kinematics 

The stellar kinematics of the eastern nucleus are more difficult to interpret; although a strong velocity 
gradient is clearly present, pinning down a rotational major axis is much less obvious. The CO(2-l) kine- 
matics are even more intriguing; the velocity map appears as a step-function, with the velocity changing 
abruptly from -135 to +205 km s^ 1 . In order to investigate this, and any differences to the western nucleus, 
we produced P-V diagrams across each nucleus, also mapping the 1.3 mm continuum emission. As can be 
seen (Fig. [T3")) . the P-V diagram of the western nucleus is fairly typical of a rotating disk or spheroid, with 
only a slight flux gap in velocity-space near the centre of rotation. The eastern nucleus however appears 
markedly different, with a gap of almost 200kms _1 between the two emission regions. To understand this, 
we construct a model of the eastern nucleus, using the simple ansatz that we are seeing two unresolved 
sources at different velocities (we note that at this stage we do not make any assumptions about the actual 
gas distribution, only the actually observed emission - as we discuss below, the gas is in fact most likely 
in a disk, but we can only see emission from two unresolved sources). After convolving a model data cube 
with the PSF of the CO observations, we extract a flux map and P-V diagram and compare it to the data. 
As Fig. [TJ] shows, the match achieved is striking - supporting our hypothesis that we are indeed seeing two 
compact, unresolved emission sources. Moreover, as Fig. 1131 shows, the peak of the continuum emission lies 
at the centre of rotation, and hence the gravitational potential. We constrain the mass enclosed, by assuming 
that the two sources are on circular orbits, and calculating the enclosed mass for a range of inclinations of 
the rotation plane. Excluding the most extreme inclinations (q<0.25 & q>0.9), we derive a dynamical 
mass of (1.8±0.5) xl0 9 M Q within 0.22", or 81 pc. However, this must be regarded as a lower limit (unless 
the gas is in a perfectly thin disk), since it neglects dispersion. 

This appears strongly suggestive of a gas disk centred on the continuum emission. However, one might 
wonder why in this case we are not seeing strong CO(2-l) emission from th e centre. There clearl y are 
significant amounts of gas present in the line of sight to the continuum peak - ISakamoto et al. I (120091) find 



the HCO + emission from the centre to display P-Cygni profiles. This, as they elaborate, is indicative of 
outflowing gas absorbing the continuum emission. We confirm the presence of absorbing gas by noting 
that the (continuum-subtracted) CO line emission at velocities near systemic appears ring-like around the 
continuum peak (Fig. I15p - a clear signature of continuum-absorption. It is therefore plausible that colder 
gas in the line of sight also causes a drop in CO line emission through line self-absorption. This also implies 
that the geometry of the eastern nucleus is similar to that of the western nucleus, in that there is a very 
dense, compact obscured central region about 50 pc across which is completely obscured from the far-IR to 
shorter wavelengths; and that around this there is a smaller, but still significant, mass distributed over a 
large region with much lower mass surface density, and hence much lower obscuration. 



6.2. Dynamical Mass Estimate 

The stellar kinematics of the eastern nucleus do not allow a straightforward measurement of a rotation 
major axis, we therefore instead use the kinematic centre and PA of the CO(2-l) velocity field (Fig. [5]) to 
measure the stellar rotational velocity and dispersion at a radius of ~ 230 pc (100 pc). Analogous to tj5.3l this 
yields an enclosed mass of 11.2 x 10 9 M Q (5.8 x 10 9 M Q ). As outlined in tj6.ll the CO kinematics suggest a 
lower limit to the dynamical mass of (1.8±0.5) x 1O 9 M within 81 pc. We note that this indicates very similar 
masses for western and eastern nucleus within 100 pc radius, and, due to the more centrally concentrated 
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distribution, that at smaller radii the e nclose d mass of the eastern nucleus is in fact larger than that of the 
western nucleus. Downes fc Solomon ( 19981 ) estimate a gas mass of 1.1 x 10 9 M Q for the eastern nucleus, 
again implying a gas fraction of ~ 10%. 



6.3. Nature of Eastern Nucleus 

As our data show, differ ent tracers display different morphologies at the eastern nucleus (cf. also the 



high-resolution Keck maps of lSoifer et al.lll999j ): raising the question whether we are seeing a kinematically 
coherent structure, or something different (a gas streamer, e.g.). The 2D velocity fields of the CO(2-l) 
and various K-band tracers all appear broadly similar, suggesting that they are all governed by the same 
gravitational potential. One way to look at this is to directly compare the velocities of the different tracers. 
Fig. [16] plots the velocity of the emission peak of a number of tracers along the CO(2-l) rotation major 
axis; as can be seen, the velocity differential is largest for CO(2-l), followed by the K-band gas tracers, and 
somewhat smaller still for the stars. Since the 2D velocity fields of these all appear qualitatively broadly 
similar, these differences are most likely due to the different depths probed by the CO(2-l) and the K-band 
tracers, rather than inherently different 3D velocity structures. The slightly smaller stellar velocities are 
expected, since rotating gas quickly collapses to a disk, whereas the stars are likely in a more spheroidal 
distribution, and hence have a larger fraction of their kinetic energy locked up in dispersion, rather than 
rotation. We therefore conclude that the eastern nucleus, like the western nucleus, is a kinematically coherent 
structure. This, in conjunction with our stellar population analysis (33])> leads to the obvious conclusion that 
most likely the nuclei are (remnants of) the progenitor galaxies' bulges. 

Since mm-wavelength CO(2-l) observations are significantly less affected by extinction than K-band 
tracers, and given that our analysis above clearly shows that our CO(2-l) observations are probing deeply 
into the eastern nucleus, we can use these to assess the depth to which our K-band tracers are probing. We 
do this by utilising our model from SjG.ll again. and making the simplifying assumption that the gas and stars 
we are observing in the K-band are also moving in a (thick) disk structure in the same gravitational potential 
as the cold gas. From this starting assumption, we then calculate at what distance above the disk midplane 
the different K-band emitters must be rotating to have velocities as observed. We find that generally they 
are 0.5-1 times as high above the disk midplane as their radial distance from the rotation axis - implying 
that we are only probing the outer regions of the eastern nucleus in the K-band, in agreement with the large 
extinction we found earlier ($4|). 

Fig. [T7] shows the extinction map, with the emission peaks of the various tracers, and the regions of 
high EWB r7 (> 16 A, black) and EWco (> 16 A, white) which indicate regions with large flux contributions 
from recent star formation. As can be seen, both the region with lots of recent star formation and the region 
with high extinction are elongated along the major axis of the CO(2-l) rotation, but offset from each other, 
with one either side of the major axis. This suggests that we are seeing a thick disk/spheroid embedded in 
the larger scale gas disk, with the northwestern half in front, and the southeastern half behind the gas disk; 
leading to the offset regions of high extinction (lower half of nucleus obscured by gas disk) and high star 
formation (nucleus above disk plane). We have indicated the sugg ested outline of the nucleus in Fig. [TTj 



We note that this is supported b y Fig. 8 of ISakamoto et al.l (|2008l ). which shows the location of the radio 
supernovae/supernova remnants (jParra et al ■I l2007t ) in relation to the 860 /im continuum emission 
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7. Nuclei &; Large-Scale Gas Disk 



As IScoville et al.l (|1997|) , iDownes fc Solomon! (|1998h , and iMundell et al.l (120011) show, the two nuclei ar e 
embedded in a larger-scale gas disk, with an estimated gas mass of ~3 x 10 9 M Q ( Downes fc Solomon 1998h . 
Our observations are consistent with this; the H2 kinematics clearly show rotational motion across the full 
SINFONI field of view (Fig.|J). 

As ou tlined in ^5]& &\ the western nucleus is quite compact, whereas the eastern nucleus appears more 



extended. IScoville et al.l (|1998l ) propose that the crescent-like K-band continuum morphology of the western 



nucleus is due to an embedded opaque dust disk, which absorbs the stellar light from the southern half of 
the nucleus (their Fig. G). This is supported by the stellar velocity map presented here (Fig. 3,j), which shows 
that the stellar rotation field is regular and exten ds to the proposedly extincted lower half. The embedded 
dust disk appears limited in extent, since in the IScoville et al.l (|1998) NICMOS data, emission from the 
southernmost tip of the western nucleus is visible. Turning to the eastern nucleus; as we argue in tj6.31 the 
extinction pattern and EWco suggest that here we are seeing a more extended spheroid, the south-eastern 
half of which is more strongly redde ned due to bein g embe dded in the larger-scale gas disk. While confirming 
the general picture put forward by IMundell et al.l (| 2 lh (cf. their Fig. 7) , in which the eastern nucleus is 
coplanar with the main CO disk, and the western nucleus lies above it, we have shown that the eastern 
nucleus is half embedded in the main CO disk. Furthermore, as our CO and stellar kinematics maps show, 
the relative velocities of the two disks are nearly identical. This, and the roughly face-on orientation of the 
large-scale gas disk, suggests that the merger is taking place roughly in the plane of the sky. 



8. Conclusions &: Summary 

We present new adaptive optics integral field spectroscopy near-infared data, and re-analyse interfero- 
metric mm CO(2-l) and continuum emission observations of the prototypical ULIRG Arp220. Our main 
conclusions are: 

• We show that in the central kpc, a 10 Myr old starburst provides the majority of the K-band luminosity, 
with a > 1 Gyr old stellar population accounting for the remainder of the luminosity and accounting 
for much of the mass. There is also a small contribution by stars < 8 Myr old, which are responsible 
for the Br7 emission. 



The ICalzetti et al.l <|2000l ) reddening law is found to provide the best fit to photometric data points 



spanning 0.45 fim to 2.12 /im. 

• Estimating the system's bolometric luminosity from IRAS fluxes, we find that we only see ~ 7% of 
the observed (non-reddening corrected) K-band luminosity. Since a reddening correction would only 
increase the observed flux by a factor ~ 2-3, this implies that there is a significant amount of emitted 
K-band light which is completely obscured. 

• Assuming the CO(2-l) emission and the K-band continuum peaks of the western nucleus to coincide, 
we find that the stellar kinematic centre is coincident with the mm continuum emission peak. This 
indicates that the mm continuum emission marks the centre of the gravitational potential in which the 
stars and gas are moving. The CO(2-l) and the stellar kinematics agree well, and all the K-band tracers 
peak at the same position; this indicates that the western nucleus is a compact, coherent structure. 
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• The eastern nucleus displays a markedly different morphology. Two emission peaks are seen in CO(2- 
1). We model the CO(2-l) data as two unresolved sources orbiting around the centre of continuum 
emission; this produces an excellent match. We find a clear signature for continuum absorption by the 
CO line-emitting gas; it is plausible that self- absorption also leads to the curious CO(2-l) morphology. 
Comparison of the velocities of the mm and K-band tracers shows that the CO(2-l) has a larger 
velocity gradient, which is expected since mm observations probe deeper into the obscured nucleus. 
We therefore conclude that the eastern nucleus also is a coherent structure, albeit more extended than 
the western nucleus. This, in conjunction with our stellar population analysis, strongly suggests that 
the nuclei are (remnants of) the progenitors' bulges. 

• This indicates that both nuclei have a similar structure in that there is a very dense and massive, but 
compact, structure in the central ~ 50 pc which is almost completely obscured at wavelengths shorter 
than far-IR; and around this is a more extended component with much lower mass surface density and 
lower extinction, but which is r esponsible for most of the observed luminosity. In general, we agree 



with the geometry suggested by iMundell et al.l ([20011 ). and can refine it further. Comparison of the 



extinction map, EWco and EWs r7 suggests that the lower half of the eastern nucleus is much more 
obscured than the upper half. This can plausibly be explained, if the eastern nucleus is coplanar with 
the larger-scale gas disk, and its lower half is buried in, and hence observed through, this gas disk. The 
very similar systemic velocities of the two nuclei, and the nearly face-on orientation of the large-scale 
gas disk suggest that the merger plane is approximately face-on. 
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Table 1. Stellar Populations 





Lk & 


EW CO a 


EW Br7 a 


LK,non-SB/LK h 


L bo i/L K h 


LK,non-SB/ Lk C 


Lboi/LK c 


SINFONI FoV 


3.25xl0 9 L© 


14.78A 


7.09A 


0.13 


151±25 


0.44 


49±9 


western nucleus 


9.5xl0 8 L o 


15.72A 


8.69A 


0.08 


111±12 


0.26 


51±5 



a measured 

b assuming the non-starburst population to be > 8Myr old with EWco=0 A 
c assuming the non-starburst population to be < 20Myr old with EWco= 12 A 

Note. — Stellar population analysis. We assume the starburst population to have EWco=17A, as measured 
between the two nuclei, and calculate the luminosity contribution of a younger (EWco=0A) or older (EWco= 12 A) 
stellar population required to dilute the CO absorption bandheads to the measured values. We furthermore calculate 
the expected L^oi/Lk for each case. The EWs r7 furthermore indicates the presence of a very young (< 7Myr) stellar 
population, the luminosity contribution of which however, at less than a few percent of the total, is negligible. 
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Fig. 1. — Spectra extracted in 0.25" radius apertures on the western nucleus (WN), eastern nucleus (EN), 
and in between the two nuclei (central). 
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Fig. 2. — We estimate the PSF of our AO data by finding the broadening function which, applied to a higher- 
resolution image with known PSF, best matches the AO image. The AO PSF is then estimated by convolving 
this broadening function with the PSF of t he high-resolution im age. Left: Our SINFONI K-band image. 
Middle: Archival NICMOS K-band image jScoville et al.lll998l ). rotated and resampled to the SINFONI 
pixel scale. Right: Best-matching broadened NICMOS image. We find the PSF to be well represented by a 
slightly asymmetric Gaussian with major axis oriented 4° east of north and FWHM 0.30" x 0.31". 




Fig. 3. — Stellar continuum, velocity, velocity dispersion, and EWco maps. White contours trace the 
continuum. The bar indicates 1", the PSF is displayed in the lower left of panel 1. North is up and east is 
to the left. 
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Fig. 4. — Flux, equivalent width, velocity, and velocity dispersion maps for a number of H- and K-band 
tracers (cf. iJ2.1.5|) . 
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Fig. 5. — C0(2-l) flux (left, with contours showing the K-band continuum), and velocity map (right, contours 
tracing the CO (2-1) flux). The bar indicates 1". North is up and east is to the left. 
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Fig. 7. — Evolution of EWco, EWsr-7, Lboi/Lx, and M/Lk, computed with STARS, for an instantaneously 
decaying starburst (red), starbursts with decay timescale lOMyr (green) and 100 Myr (blue), and continuous 
star formation (black). As can be seen, only a ~ 10 Myr old 'delta'-starburst can produce EWco ^ 14 A. 




Fig. 8. — Comparison of different extinction models: We fitted a synthesised stellar population spectrum to 
HST photometric data points extracted in 0.15" radius apertures in a range of locations across our field of 
view (4 examples shown), by redde ning the synthetic SE P until it matched the data. Grey: mi xed extinction, 
black: screen extinction, magenta: Calzetti et al. ( 2000l ) reddening law. As can be seen, the ICalzetti et all 
((2000) reddening law produces the best fit to the data. 
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Fig. 9. — Left: F b s /F em in the K-band, derived by fitting a stellar template to the line-free continuum using 
the ICalzetti et al.l |200d ) reddening law. Right: corresponding optical extinction Ay. Contours trace the 
continuum. 




Fig. 10. — Left: F b s /F em in the K-band, derive d by comparing the o bserved flux ratio of Paa and Br7 to 
the theoretically expected value, again using the Calzetti et al. ( 2000h reddening law. Right: corresponding 
optical extinction Ay. Contours trace the continuum. 




Fig. 11. — Western Nucleus: Stellar velocity field, best- fitting disk model, and residuals. The kinematic 
centre is marked by a black circle, and the major axis is indicated with a white dashed line. 1 pixel corresponds 
to 0.05". 
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Fig. 12. — Western Nucleus: C0(2-l) velocity field, best-fitting disk model, and residuals. The kinematic 
centre is marked by a black circle, and the major axis is indicated with a white dashed line. 1 pixel corresponds 
to 0.05". 
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Fig. 13. — Position- Velocity diagrams across the western and eastern nucleus. On the left is a CO(2-l) 
flux map indicating the position of the slits, with 1.3 mm continuum emission overlaid in contours. On the 
right are the P-V diagrams for western (top) and eastern (bottom) nucleus. As can be seen, there is no line 
emission in the central region of the eastern nucleus, coincident with the peak in dust emission. 
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Fig. 14. — Position- Velocity diagram across the eastern nucleus (solid black contours) and from model 
(dashed red contours). Both sets of contours are displaying the same flux levels. As can be seen, the model 
(two unresolved sources orbiting the continuum emission peak) provides an excellent fit to the data. 




Fig. 15. — Integrated CO(2-l) flux map, with insets showing the CO line flux from a narrow velocity 
range near systemic velocity; contours mark mm-continuum emission. As can be seen, in those wavelength 
intervals, the (continuum-subtracted) line emission drops markedly towards the continuum peaks. This is a 
clear indication for continuum absorption by the line-emitting gas. 
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Fig. 16. — Position and velocities of CO(2-l) and various K-band tracers in the eastern nucleus, along the 
CO(2-l) rotation major axis. CO(2-l) displays the largest velocity differential, followed by the K-band gas 
tracers, and finally the stars. This is interpreted as arising due to the different depths probed by the various 
tracers. 
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Fig. 17. — Extinction map (cf. Fig. [9J, with the emission peaks of the various tracers, and the regions of 
high EWfi r7 (> 16 A, black contours) and EWco (> 16 A, white contours). The latter indicate regions with 
large flux contributions from recent star formation. As can be seen, both the region with lots of recent star 
formation and the region with high extinction are elongated along the major axis of the CO(2-l) rotation, 
but offset from each other, with one either side of the major axis. This suggests that we are seeing a thick 
disk/spheroid embedded in the larger scale gas disk, with the northwestern half above, and the southeastern 
half below the gas disk; leading to the offset regions of high extinction (nucleus obscured by gas disk) and 
high star formation (nucleus above disk plane). The black dashed line is a suggested outline of the nucleus. 



